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Page 7; In equation (13) the plus sign preceding the last term (the 
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This last item wenrants the following explanation: 
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was ottained by evaluating the integreil 
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by con 5 >arison with the known value of the integral 



This procedure leads to an incorrect result since the branch points 
of the Hankel functions involved are on the real axis in one case nnt^ on 
the imaginary axis in the other. 

The integral being considered actually appears in an infinite 
summation of integrals, namely 



A correct value of this summation is obtained by replacing the series 
of Hankel functions in the last integral by the equivalent series of 
exponential functions given in eqiiation (l8) of the report and th^n 
performing the indicated integrations. 
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This procedxire leads to the following final expression for 

(eq. (17)) 
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CONSIDERATIONS ON THE EFFECT OF WIND-TUNNEL WALLS ON 
OSCILLATING AIR FORCES FOR TWO-DIMENSIONAL 
SUBSONIC COMPRESSIBLE FLOW 
By Harry L. Rvinyan and Charles E. Watkins 


SUMMARY 


This paper treats the effect of wind-tunnel walls on the oscillating 
two-dimensional air forces in a coii?)ressible medium. The walls are simu- 
lated by the usual method of placing images at appropriate distances 
above a,nd below the wing. An important result, shown is that, for pertain 
conditions of wing frequency, tunnel height, and Mach number, the tunnel 
and wing may form a resonant system so that the forces on the wing .are 
greatly changed from the condition of no tuiuiel walls,.. It is pointed 
out that similar conditions exist for three-dimensional flow in circular' 
and rectangTilar tvinnels and apparently, within certain. Mach number ranges, 
in tionnels of nonuniform cross section or even in open .tunnels or jets. 


iNTRODUCnON 


The understanding of flutter and other nonsteady phenomena requires 
a knowledge of the associated unsteady flow. In the underlyuig theories 
of unsteady flow, such assximptions as small displacements, linearizations, 
and an inviscid fluid are made in order to obtain workable and usable 
resiilts. When it is necessary to investigate the effect of these assump- 
•^ions on analytical results by measurements of the forces and moments on 
an oscillating wing in a wind tunnel or to treat cases that do not con- 
form to theory, the question of the effect of the tunnel walls naturally 
arises. In the case of steady flow the problem of the effect of tunnel 
walls is more or less classic and has been treated by many investigators. 
In general, these investigators have been able to obtain relatively simple 
factors which can be used to modify measiirements of the air forces on a 
wing in a tvinnel to correspond to free-air cpnditions. The’ extension . 
of the results to compressible flow presents no difficulties since the 
results for incoipressible flow can be corrected according to Prandtl— 
Glauert correction factors . 
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In the case of unsteady flow, Reissner, reference 1 , and W. P. Jones 
reference 2 , have published papers showing the effect of wind-tunnel walls 
for the incompressible .case. In both papers, the influence of the tunnel 
walls is foimd to be comparatively small for most cases, although indica- 
tions are given that, for some ranges of a reduced-frequency parameter, 
the effect may be quite large. In the unsteady case, \inlike the steady 
case, the transition from resiolts for incompressible flow to those for 
coirpressible flow cannot be accomplished by sin^^le transformations. This 
difficulty is a result of the fact that, in an incompressible fluid, the 
velocity of propagation of a disturbance is infinite and no time lag 
occi^s between the initiation of a disturbance and its effect at another 
position in the field, but, in a compressible flxiid, a definite time is 
required for a signal to reach a distant field point so that both a phase 
lag and a change in magnitvide result. Under certain conditions this 
phase lag can result in a resonant condition which would involve large 
corrections . 

The puip)ose of this paper is to consider the effect of wind-tunnel 
walls on the forces on an oscillating airfoil of infinite span with con- 
siderations of the compressibility of the flviid. The usual method of 
inages is employed in order to satisfy the condition of no normal velocity 
at the tunnel walls. First, the effect of tunnel walls on the induced 
vertical velocity, hereinafter referred to as downwash, of an oscillating 
doublet is determined and this result is used to formulate the integral 
equation for the downwash of an oscillating airfoil in a tunnel. This 
paper is npt intended to give numerical values or any detailed calcula- 

of final tunnel-wall correction factors but mainly to show the 
existing need for such calcxxlations and to present equations for calcu- 
lating corrections for the two-dimensional case. 


SYMBOLS 


A 


constant 


b semichord 

Q velocity of sound 

H tuinnel height 

Hankel fimctions 
M Mach number 

4 p local pressure difference 
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V 


w, Wq, 

I, X, y, z, 
P = /l - 
7 

CO 


X 

9 

P 


time 


velocity 

dovmwash or vertical induced velocity 
Cartesian coordinates 

Euler's constant 
angular frequency 
wave length 
acceleration potential 
velocity potential 
fluid density 


ANALYSIS 

Effect of Tunnel Walls on the Downwash of a Single Doublet 


The differential equation that governs flow due tc) small ndhsteady 
perturbations imposed on a steady, uniform flow field is the wave equa- 
tion. Referred to rectangular coordinates, fixed relative to the vindis- 
turbed stream at infinityj this equation is 


^ ^ ^ _ i. ^ = 

3y^ ® c^ 3t^ 


( 1 ) 


In this equation the independent variable ijr may be regarded as either 
a perturbation velocity potential or as an acceleration potential* In 
treating the boundary conditions of the second section of this analysis 
it is convenient to regard ijr as an acceleration potential. Thus, in 
order to be consistent, i is hereinafter regarded as an acceleration 
potential. Accordingly ilr is directly proportional to a pertxirbation 
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pressure field and is therefore related to a perturbation velocity 
potential cp as follows: 


i = 


3cp 3cp 

— + V — 

3t 3x 


( 2 ) 


In order to calculate the downwash w = associated with it is 

dy* 

necessary to solve equation (2) for 9 in terms of 

When \jr and 9 are sinusoidal functions of time, such that 


'Kx,y,t) 

9(x,y,t) 



( 3 ) 


equation (2) becomes independent of time and thus reduces to an equation 
with one dependent variable, namely 




IC09 


V ^ 
3x 


This equation can be integrated id.th respect to x to .give 


Oi) 



( 5 ) 


where the lower limit of integration is chosen for later convenience so 
that 9 vanishes far ahead of the point of disturbance. The downwash 
may be readily calculated with. the use of this equation. In the absence 
of tunnel walls the retarded potential (that is, the potential corr 

spondi^ to outgoing waves) of a harmonically pulsating pressure doublet 
located, for simplicity, at (0,0) that satisfies equation (l) is 
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where and are Hankel functions as defined in reference 3, 

A is an arbitrary c onstant denoting doublet strength, co is circular 

frequency, and p = j/l - The Hankel fimction in equation (6) 

becomes infinite ^s ^ ^ j as its argument approaches zero. Other- 

( 2 ) . ^ 

wise is continuous and approaches • zero as its argument approaches 

infinity. Thus the only discontinuity in \|tq is at the location of the 
doublet, that is, at (x=0,y=0). 

In the presence of plane tunnel walls located parallel to the x-axis 
at H/2 units above and H/2 omits below the doublet position, the 
potential \|^ of a pressoire doublet may be represented by the potential 
of an infinite system of appropriately chosen reflecting doublets, namely 
(see fig. 1) 


= Airie 




r 


cp^ * 


n— CO 


(7) 


In this equation the term corresponding to n = 0 is the potential fg, 

equation (6), discussed in the preceding paragraph. It may be noted that 
only this temi of the infinite summation in equation (7) gives rise to a 

discontinuity in f at any point within the tiinnel -oo <X<ce\ 
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The . inf ini ty of terms corresponding to n 0 are necessary to cause 
the downwash w to vanish at all points of the tunnel walls, y = ±S. 

The downwash along the midsection of the tiinnel y = 0 is given by 



represents the downwash associated with the pressure doublet in the 
absence of tunnel walls and 



( 10 ) 

represents the additional downwash due to the presence of tunnel walls. 
Thus the relative value of wg as con^ared with wg + w^ is the main 
item of interest here. 

The integrals appearing in equations (9) and (10) can be reduced to 
sin^jler form for evaluation but since the steps required to reduce one 
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of the integrals are the same as required to reduce the other, only 
the integral appearing in equation (9) will be treated in detail. The 
reduced form of the other integral can then be obtained by simple com- 
parison, The Hankel function in equation (9) satisfies the following 
identity: 



Substituting this relation into equation (9) gives 



In equation (l?) the first integral can be . integrated twice by parts to 
give for wq ^ 
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By vjriting the integral in equation (I3) as the sum of two integrals, 
namely 



and making a change of variable 


col 

Vp2 


= u 


the expression for Wq may be further reduced to 


Wrv = 


Am 


lim e 
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‘iMe^P‘ 


y->0 
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(15) 


( 16 ) 




In the limit y = 0 the expression in braces in equation (16) reduces 
to the kernel of Possio's integral equation relating pressure and down- 
wash for the, oscillating airfoil in cor^jressible flow. (This result 
checks, the results for this expression given, for exait^Dle, in reference U.) 



The value of the integrals in equation (10) may be similarily reduced to give 
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In general, this infinite-series representation of Wq_, equation (I7), 

converges to a finite value. However, for certain critical values of the 
frequency parameter coH/V, it is found that the value of becomes 

infinite. This fact can be readily made evident by use of relations 
given in reference 5 where it is shown that an infinite series of Hankel 
fimctions of the type appearing in equation (17) can be replaced by an 
equivalent series of exponential functions as follows: 



It may be seen that, if this relation is substituted into equation (I7), 
the value' of W]^ becomes infinite for all values of x when the fre- 
quency parameter coH/V has any of the values given by 

f - ( 2 m - 1) ^ (m = 1 , 2 , 3, . . .) ( 19 ) 
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These critical values of the frequency parameter correspond to a condi- 
tion of pure resonance in the tunnel which in the present case implies 
that a harmonic disturbance of any finite amplitude may lead to a down- 
wash of infinite anqjlitude. 

Of co\irse these infinite values of w^ would never be realized 
\inder practicable conditions because factors such as finite tunnel length, 
absorption through walls, fluid viscosity, and so forth that would give 
rise to dan^jing would make pure resonance unobtainable; however, with 
damping present, resonant frequencies yielding values of ccH/V in the 
neighborhood of those given in equation (19) would exist and it is not 
likely that quantitative agreement or even possibly qualitative agreement 
between calculated and measured dovmwash (or forces) can be realized when 
the value of coH/V is in the neighborhood of these critical values. 

It is interesting to note that the effect of boundary conditions 
such as section geometry, tunnel-wall flexibility, and so forth is to 
change the value of the critical frequency but not to do away with the 
possibility of resonance. Also, by treatments similar to those eii 5 )loyed 
herein, it can be shown that under idealized conditions resonance can 
occur in three-dimensional flow in both round and rectangular tunnels or 
apparently, within certain Mach number ranges, in tunnels of nonuniform 
cross section (expanding or contracting section) or even in open tunnels 
or jets.. 

The fundamental or smallest critical values of ccH/V, corresponding 
to m = 1 in equation (19), are shown plotted as functions of Mach num- 
ber M in figvire 2. This figure indicates that there is no finite 
critical value of ocH/V for the conditions M = 0, V ^ 0, and c = », 
which' correspond to a flow of incompressible fluid in the turinel. This 
result agrees with those found in references 1 and 2 . 

The frequency parameter 


SS = (2m - l)np (m =1, 2, 3, . , .) (20) 

c 

which may be derived from equation (19) .is shown plotted, for m = 1, 
as a function of Mach number in figure 3. Equations (20) and figiire 3 
show that finite values of the critical frequency exist for the condi- 
tions M = 0, V = 0, and c ». These conditions correspond to a com- 
pressible fluid at zero velocity in the tunnel. For these conditions 
equations (20) and the corresponding wave lengths 



2H 



2m - 1 


(m 1, 2, 3 , • 


( 21 ) 
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agree, respectively, with results found in the literature for the char- 
acteristic frequencies and wave lengths associated with transverse 
acoustic vibrations in rectangular chambers when the location of the 
source of disturbance is excluded as a nodal point. See, for exaiT 5 )le, 
reference 6. 

It may be of interest to note that equation (20) can be derived 
from the principle of standing waves as follows: The condition for 

resonance for the type of distxirbance considered in 5 )lies that the 
standing transverse waves have a maximum velocity at the midsection of 
the tunnel and zero velocity at the boundaries. A half -sine wave of 
wave length X = 2H or any odd divisor of this length, namely, 

2H 

X = satisfies this condition. If c is the velocity of sound 

2m - 1 

in the medium and V the velocity of the medium, the velocity of pro- 
paga tion of a distxirbance in a fixed plane perpendicular to the air flow 

is /c^ - V^. Since the frequency is given by the speed Of propagation 
divided by the wave length there is obtained 

f = - l)/c^ - ^ 

2H 2tt 


or 


f = tfl(2m - 1) 


Integral Equation for an Airfoil of Infinite Aspect 
Ratio Oscillating in a Wind Tunnel 

In order to present equations from which tunnel-wall corrections 
for two-dimensional flow can be calculated, use is made of the foregoing 
analysis to derive the integral equation, relating downwash distributions 
and lift distributions, for the effect of tunnel walls on the lift dis- 
tribution associated with a given downwash distribution. 

The resultant pressure or local lift Ap associated with the accele- 
ration potential of a single doublet located at (xq,0) with strength 

depending on streamwise position Xq may be expressed simply as (compare 
with equation (6)): 
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For an airfoil inside a two-dimensional tunnel the relation between downwash and local lift 
becomes (compare with equation (8)) 
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It appears desirable to solve equations (23) and (2I4.) by collocation or some other approxi 
mate method to obtain tunnel-wall corrections for some particular cases of prescribed downwash 
and to determine experimentally the range, if any> of frequency parameter in which quantitative 
resiilts can be obtained for these cases. 
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CONCLUDING REMARKS ' 


The important result shown is that, in a tunnel of infinite length 
containing a flowing fluid, a resonant condition involving a transverse 
oscillation of the fliiid across the tunnel is possible and measured air 
forces at or near this condition of resonance might be greatly modified 
from those measured in free air. This resonant condition is a (sin^jle) 
function of Mach number, timnel height, and wing frequency and brings to 
attention a new type of tunnel-wall interference. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va., September 2li, 1951 
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Figure 1 .- Sketch showing reflecting system of doublets sim\ilatlng 

two-dimensional txuinel walls. 
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Figure 3 .- Fundamental critical values of frequency parameter oil/c 
plotted as a function of Mach number M, 
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